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Abstract 

We consider large-order perturbative expansions in QED and QCD. The coefficients of the 
expansions are known to be dominated by the so called ultraviolet (UV) renormalons which 
arise from inserting a chain of vacuum-polarization graphs into photonic (gluonic) lines. In 
large orders the contribution is associated with virtual momenta A;^ of order Q^e" where 
Q is external momentum, e is the base of natural logs and n is the order of perturbation 
theory considered. To evaluate the UV renormalon we develop formalism of operator product 
expansion (OPE) which utilizes the observation that ^ Q^. When applied to the simplest 
graphs the formalism reproduces the known results in a compact form. In more generality, 
the formalism reveals the fact that the class of the renormalon-type graphs is not well defined. 
In particular, graphs with extra vacuum-polarization chains are not suppressed. The reason 
is that while inclusion of extra chains lowers the power of In fc^ their contribution is enhanced 
by combinatorial factors. 
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1 Introduction 



Large-order behavior of perturbative expansions has been studied for about forty years start- 
ing from the the seminal paper by Dyson (for review and further references see Refs. 

|]). Ahhough a lot of insight has been gained some practical aspects of the issue in 
case of QCD have not been clarified so far and attracted for this reason a renewed attention 
recently. In particular, let us mention here discussion of possible 1/s power corrections at 
large energies a/s, as generated through divergences of perturbative expansions [§, H, 0- 
In this note Q we address the problems of calculation and of calculability of the so called 
ultraviolet renormalon p, ^ |10|, [Tl| which are related to these 1/s corrections but of more 
general nature. 

The importance of the ultraviolet renormalon rests on the observation that it dominates 
large orders both in QCD and QED To be more specific we shall have in mind pertur- 
bative calculations of the famous ratio R : 

, (j(e"'"e" — > hadrons) 

Ri^) = 1 



as function of energy and its Euclidean counter-part Il{Q 



2^ 



^ dQ^ 127r2 Jo (s + Q2^2 ^ > 



as function of Euclidean momentum Q. The quantity Q^dIl{Q'^)/dQ^ is represented as an 
expansion in QCD coupling constant ag^Q"^) : 

- Q'—^ = {Pcirton model) x ^ a^a^iQ^) , (3) 

dQ n=0 



where first three coefficients have been calculated explicitly JT2|. 

Note that the normalization of the coupling at , as above, guarantees that there are 
no log factors in the expansion (|]). Since we consider Q"^ to be large the coupling ^^(Q^) is 
small. For the sake of simplicity we will illustrate our approach mostly on example of U{1) 
model where strong interactions are mediated by an abelian vector field. Then we consider 
energy to be far below the position of the Landau pole so that ai is small. As we shall see 
the conclusions are of general nature and apply both to the U{1) model and QCD. 

Now, the generic behavior of the expansion coefficients at large n looks as 

an ™ Kn^^^ (4) 



^The letter-type version of the paper was published in the Physical Review Letters, 73, 1207 (1994) 
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where K, 7 and S are constants and actually there exists a variety of sources for the factorial 
growth resulting in different S. The constant S takes the smallest absolute value in case 
of the ultraviolet renormalon: 

SuV renorm ~i (5) 

where ho is the first coefficient in the corresponding (3 function: 

= -hoa\Q') - ha\Q') + ... (6) 

Note that in QCD is positive and the series (§) is sign alternating while in U{1) case 
is negative. The sign alternation allows for the Borel summation of the series and from this 
point of view there exists an important difference between QCD and U{1). In this paper 
we are concerned mostly with calculating of a„ at large n, not on the summation of the 
perturbative series. 

The main technical tool we are using is the operator product expansion (OPE), or the 
expansion in inverse powers of k"^ where k is the virtual momentum flowing through the 



gluonic line (see Fig. 1). The idea of using this OPE is outlined first by Parisi ||T0[ and 
is based on the observation that effectively k"^ ^ in case of the ultraviolet renormalon. 
More specifically, we shall see later that 

{k\}f - Q'e\ (7) 

where n is the order of perturbative expansion and e is the base of natural logs. We will 
elaborate the idea of the expansion in jk^ on explicit examples. 

Our primary objective in this paper is to work out a scheme for explicit evaluation of 
renormalon contributions. In Sec. 2 we will address computation of the simplest renormalon- 
type graphs represented on Fig. 1. In f/(l) case such calculations have been performed in a 
number of papers |T3|, H, [1^ . In particular we reproduce the results of Beneke who 



evaluated the contribution of the renormalon chain directly, without use of OPE. What we 
would like to add here is a simplification of the scheme. Moreover, the generalization to the 
QCD case is imminent. As the next step we generalize the procedure to the case of two 
renormalon chains, i.e. three-loop skeleton graphs (Sec. 3). We will demonstrate that this 
three-loop contribution dominates over the two-loop one. 

The consideration of graphs with one or two renormalon chains naturally brings us to the 
problem of calculability of the renormalon contribution in general (Sec. 4). By this we mean 
the problem of identifying the graphs which control the coefficient in front of the factorial 
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(see Eq. (|D). In particular, one may increase the number of renormalon chains and the 
question is whether this contribution is suppressed or not. f\ 

Superficially, the more complicated graphs can be neglected since at least one extra 
factor of a is not accompanied in this case by a log. The use of the operator product 
expansion mentioned above allows to analyze the problem in a general way. The result turns 
unexpected: the contribution of the graphs with extra chains is not suppressed and the class 
of the "renormalon-type" graphs is not well defined at all. 

The reason is that we have in fact two large parameters, namely. In k'^/ and the order of 
the perturbative expansion, n. Moreover, effectively Infc^/Q^ ~ n, as is implied by Eq. (|3T|) . 
While the coefficient in front of the highest power of the log can be found in a straightforward 
way - and this is the essence of the renormalization group analysis - it turns to be not enough 
to evaluate the asymptotic of a^. The other contributions lose the log factors but gain extra 
factors of n because of combinatorics. 

2 Operator product expansion. Two-loop example. 

Let us consider polarization operator H/^u, 

U,, = zj dxe''^^0\T{j,ix)UOm = {q.qu - g.^q^MQ'), = -q^ (8) 
of electromagnetic current in the simplified model with Nf fermionic fields 

where Qq are the corresponding electric charges. Strong interactions are mediated by U{1) 
gluonic field 5^ and the Lagrangian of the model is 

L= -\B,,B^^^ + J2qYtd,q + 9B^jl + eA^j, (10) 

q 

where g and e are the strong and electromagnetic couplings, respectively, and 

= (11) 

Q 

^ In case of QCD because of the complexity of the gauge-fixing problems even the set of the graphs 
delegated to each type is not so well defined a priori. We shall see, however, that the technique developed 
allows to circumvent this problem as well. 
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The sum of electric charges of fermions is taken to be zero to avoid a mixing between fields 
and due to fermionic loops, 

= 0. (12) 

Q 

The simplest renormalon-type graphs are depicted in Fig. 1 and we will explain the basic 
features of the technique on this example. The sum of these graphs can be cast into the 
following form: 

e^n,.(g)ef,,er,, = -U ^'-^{l*\T\Y) (13) 



- 2 J (27r)4 P 

where we have made the Euclidean rotation in the integration over a gluon momentum 
and have introduced polarization vectors e'^^ 2) initial and final virtual photons. The 
running coupling g'^{k'^) sums up vacuum bubble insertions and has the form 

9\k') = 3'm(l + bo'-^ln^Y' (14) 



with the first f3 function coefficient 



Nf 

"» = - *F • ('^) 



The matrix element (7*|T|7*) represents the forward amplitude of gluon-photon scattering 
and the operator T is 

T = Jdxe"'^T{jl{x),j'^{0)}. (16) 

The relation ( [T3| ) is evident for the graphs of Fig. 1. We will discuss in Sec. 4 how it is 
modified for higher loops. 

By assumption - to be checked a posteriori - the momentum k fiowing through the gluon 
line is much larger than external momentum Q. Then it is logical to start by expanding in 
inverse powers of k. Thus we come to consider OPE for T-product of two gluonic currents 

T = I dx e''^ T{jlix)f\0)} = Y.^iik)Om, (17) 

where Oi are local operators. We will discuss this operator expansion in more detail in the 
next section. Here we find coefficients Cj in the tree approximation using the Schwinger 
background field technique [|1^] (for a review, see Ref. |1^). In this technique, the lines in 
the graph of Fig. 2 are understood as propagators in external electromagnetic and gluonic 
fields and the operator T takes the form: 

T = - dx{x\qiX)r^lMXm + {k^ -k) . (18) 
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Note that within the framework considered operators of the coordinate and of momentum 
Pfj_ are introduced. Moreover 

[X,,p.] = -ig^^u, [X^,X,] = 0, [p^,p.] = (19) 

The operator is defined as 

P^. = P^. + gB^ + eA^Qq (20) 
and the matrix element (^8|) is taken over eigenstates of operator X^, 

X^|x) = Xf,\x). (21) 



The construction of the OPE reduces to an expansion of Eq. ([18| ) in powers of . The 



first and second terms of the expansion vanish upon averaging over directions of 4-vector 
kfj_ and use of equations of motion, ipq = 0. Thus, the expansion starts from operators of 
dimension six, 

T = -^,T.inP'',[Pu,P,]]q + 0{k-% (22) 
Substituting the commutator 

[P^,P,] = igG^, + ieF^,Qq (23) 

we get 

T = ^ [ed^F,, Y: Q,qYq + gD^G,, ^7^9 ) + 0{k-'). (24) 

The next step is to evaluate the matrix element (7*|T|7*). The part of T containing 
D'^Gy^ will contribute only on three-loop level and will be considered in the next section. 
As for the part of T containing d'^Fy^ it immediately factorizes into 

(7lT|7*) = -1 [(^*|e5^F.^|0)(0|/|7*) + (7I/IO) (0|e9-F.,|7*) ] (25) 

where is the electromagnetic current (see Eq. (^). The matrix element of d^F^^ is trivial: 

(Olea'^F.^b*) = -e (g^eW - g^(ge(i))) . (26) 

The matrix element (7*|j/i|0) is given by the well-known one-loop graph (see Fig. 3) and is 
equal to 

(7 IJmIO) = ey-qM^^ n) J j^^^, = i^^l^g^-^^ e|,^-g^(ge^ 0) 

(27) 
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where {Q'l) is the averaged square of electric charge, 



Note that the integral over the fermionic loop has been evaluated with logarithmic accuracy. 
The upper limit of integration, ~ A;^, is implied by our OPE construction while the lower 
bound, ~ Q^, arises from account in the integrand for the external momentum Q. 
Substituting the result (ESf) for (7*|T|7*) into ([T3|) we come to 



^iQ") = - ^7?^Q" . .Tt1^7^-^?'(^') (29) 



Here = —q^ and the integration over Euclidean fc^ runs over k"^ > Q"^. We are interested 
in the expansion of n((5^) in aiiQ"^) = g"^ {Q"^) / An . Expanding Eq. ([T^) we get 

oo >2 

g\k') = A7:a,{Q')Y.(rh,a,{Q')Y\n^— (30) 

n=0 ^ 

where 6o is the first coefficient in the (3 function. Performing the integration over k"^ in the 
right-hand side of Eq. (|29| ) , 

poo rlh^ h"^ 

Q' L = n!, (31) 



we arrive at the final expression for the UV renormalon contribution: 
Differentiating this expression we get 



n(Q^) = ^|^E(-^oai(Q^)) n! . (32) 



= ^2^l-3^JS^"'(-'°"^) • ^''^ 

In other words, at large n the coefficients in ai expansion (for the definition see Eq. (|]) ) 
are given by 

an™-TrV(-M"-^! • (34) 

This result coincides with explicit calculations of Ref . P] . 

To summarize, in this section we utilized the operator product expansion to evaluate 
the graphs in Fig. 1 and compared the results with direct calculations. The advantage 
of the operator expansion is not only the compactness of the calculation but also ensuring 
automatically the gauge invariance so that the generalization to QCD case is straightforward 
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and reduces to the change in Bq. Indeed since the operator product expansion is based on a 
set of gauge invariant operators the only dependence on Infc^/Q^ arises through the use of 
Eq. ( pop (or its two-loop generalization) and we do not need even to specify the gauge fixing 
or the class of graphs involved explicitly. This technical point is especially important in case 
of non-abelian gauge theories. 



3 Comments on operator product expansion. 

In this section we combine a few simple remarks on the use of OPE exemplified in the 
preceding section. 

First, let us demonstrate that Eq. (|T3|) which relates the polarization operator II^i, and 
the forward amplitude of gluon-photon scattering is valid in any order of perturbation theory. 
The functional dependence on the gluon propagator Dafsi^x) is given by the following textbook 
formula: 



expli J (fxCint i ^ j i exp^- j d'^xd'^yrio,{x)Do,p{x - y)rip{y)^ 



(35) 

r)^(x)=0 



6ri^{x 

where Lint is the interaction Lagrangian in which the vector field A^{x) is substituted by the 
functional derivative 5 /5rj^{x) over the source rj^{x). The source rj^{x) is substituted by zero 
after the differentiation. To separate out the "hard" part of the gluon propagator D^p{x) 
which corresponds to the exchange by momentum k let us present Daf3{x) as 

D^p{x) = D^^p{x) + Dlf{x). (36) 

Then expanding Eq. ( ^5] ) to the fist order in D^p{x) we are proving the relation (|TB|) in higher 
orders. What is crucial for the proof it is the ordering of momenta with the momentum k 
as the highest one. 

The first comment is to note that our calculation justifies the assumption about dom- 
inance of large k"^. Indeed, the integrals (0) over fc^ are saturated at large n by a saddle 
point at kljj = Q^g" (see Eq. (g)). The width of the range of k"^ which contribute to the 
saddle point integration is 

^-p^ - ^ (37) 

It means, in particular, that the precise value of the lower limit in integrals (^Of ), (|3TD is of 
no importance at large n. 

Thus, expansion in Q"^ /k'^ is fully justified. Moreover, it is clear from the calculations 
above that it is the dimension of operators Oi (see Eq. (p^Tj)) which is most important. 
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Namely, for an operator of dimension d the contribution to an is proportional to 

^'''wfc^^^V ^ ((d-2)/2r+^- ^^^^ 

In the preceding section we dealt with operators of dimension ci = 6. Operators of higher 
dimensions give rise contributions which are suppressed by powers of (d — 4)/2 and can be 
safely neglected. 

However, the question naturally arises on the role of operators of dimension four. Al- 
though such operators did not appear in the tree approximation of the preceding section 
they show up via loop corrections. It is worth emphasizing therefore that these operators 
are not relevant to asymptotic of a„. 

Indeed, matrix elements of d = 4 operators over virtual photons have the form 

{l*\Ot'\l*) = A,e2(g2(e(i)e(2))-(ge«)(ge(2))) (39) 
where Ai are dimensionless and the corresponding contributions to n((5^) are 

~ j deAici=\e)g^{e) (40) 

where cf^^ are OPE coefficients, cf^'^ ~ To get rid of the apparent ultraviolet diver- 

gence in (^OD one needs to consider instead of n((5^) the quantity ■ dI{{Q'^) / dQ'^ which 
contains no ultraviolet uncertainty. 

The finiteness of ■ dIl{Q^)/dQ'^ implies that J2icf^^Ai is independent of Q^. Notice 
that this assertion is valid in spite of nonzero anomalous dimensions of some d = A operators. 
It also persists in the presence of fermionic mass terms m^gg which contains d = 3 operators. 
The independence J2i cf^^Ai in any order of perturbation theory means that the sum of 
contributions (|0]) drops off from ■ dIl{Q'^) / dQ"^ . 

Thus d = 4 operators do not contribute to UV renormalons and the UV renormalon 
calculus deals with d = 6 operators. In the preceding section we had one such operator 
explicitly: 

Ofj = ed^F.^.f (41) 

where is defined by Eq. (H). The calculation of matrix element {l*\OFj\l*) (see Eq. (p5|)) 
can be interpreted as a mixing of Ofj with operator 

Of2 = {ed-'F^^f, (42) 

i.e. with logarithmic accuracy we have 



127r2 g2 



OFjlk^ = OfjIq^ + ——^\n— ■ Of'^\q2 (43) 
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where we marked the normahzation points of the operators. Since the diagonal anomalous 
dimensions of the operators 0Fj,0p2 are vanishing Eq. ( ^3] ) specifies the matrix of the 
anomalous dimensions completely. 

Note also that the log in the right-hand side of Eq. (^) results in an enhancement of 
a„ by factor of n (see Eq. (PT])) and this enhancement is crucial for the consistency of our 
calculation. The point is that in one- loop order the operator Of2 appears in the OPE ( |T7|) 
not only through its mixing with Opj as in the tree approximation but directly as well. In 
that case large momentum k flows through all fermionic lines in graph of Fig. 1. However 
then there is no log factor similar to the one in Eq. (^) so that the direct Op2 coefficient 
can be neglected with 1/n accuracy. There are also d = 6 four-fermionic operators which 
contribute to n(Q^) on three-loop level and these will be discussed in the next section. Here 
we would like only to note that their logarithmic mixing with the operators Opj and Of2 
results in an enhancement of the three-loop skeleton graph. 

Finally we note that so far we considered for simplicity one-loop f3 function for ai. 
Accounting for the second coefficient in the j3 function brings a non-vanishing constant 7 in 
Eq. (^). The procedure is rather standard and we will give examples in the next section. 
As for the overall normalization (constant K in Eq. (^) it is actually scheme dependent, 
for further discussion see |jl8|, |19|. The constant K can be fixed in large Nj limit, which 
simplifies calculations greatly (see, e.g., [|, 0). As we shall see later, however, the limits 
of n and Nf tending to infinity are not necessarily consistent with each other. Here we keep 
track of large n dependence. 

To summarize, the operator product expansion provides with systematic means of calcu- 
lating asymptotic of the perturbative expansion coefficients a„ associated with renormalon- 
type graphs. 



4 Three-loop example. 

In this section we evaluate the renormalon contribution associated with the graphs with 
two chains of vacuum-polarization insertions, or three-loop skeleton graphs, see examples 
in Fig. 4. In fact we have already mentioned that the contribution of the operator D'^G^^ 
(see Eq. (0)) arises only on a three- loop level and now we come to consider three- loop 
calculations in more detail. 

Note that D^G^^ can be substituted by (— (?) Y.'iln'i via equations of motion so that we 
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get for the corresponding part Ti of operator T defined by Eq. (|T6D 



Ti = -%0,, Oi= (^Em) • (44) 

Diagramatically the appearance of T\ in three-loop graphs for the polarization operator is 
illustrated by the diagram h in Fig. 4. The dotted box in this diagram presents a "penguin" 
mechanism for T\ similar to the one in weak interactions . 

A different mechanism leading to four-fermion operators is due to the exchange by two 
gluons, see the diagram a in Fig. 4. The part of the graph which corresponds to the operator 
T is marked by the dotted box in Fig. 4a and presented separately in Fig. 5. As far as explicit 
calculations of the graphs are concerned they are very similar to those used in derivation of 
QCD sum rules |^ since OPE is exploited in the both cases. The difference is that in the 
former case it is the momentum carried by the electromagnetic current which is considered 
to be large while now the large momentum is brought in by the gluonic current. 

The subtle point here is that the relation between the scattering diagram of Fig. 5 and 
the box part of diagram a in Fig. 4 contains a combinatorial coefficient 1/2 which is absent 
for the "penguin" box of Fig. 46. Indeed, the functional dependence on the gluon propagator 
Dapix) is given by the following textbook formula: 

6_ 

where Lint is the interaction Lagrangian in which the vector field A^{x) is substituted by the 
functional derivative 5 /5rj^{x) over the source rj^{x). The source rj^{x) is substituted by zero 
after the differentiation. To separate out the "hard" part of the gluon propagator D^p{x) 
which corresponds to the exchange by momentum k let us present Dap{x) as 

DM = Dl,{x) + D-f\x). (46) 

Then expanding Eq. (|^) to the fist order in D^p{x) we get the relation (pIS)). To maintain 
the same relation (|13]) for terms of the second order in D^p{x) the corresponding part of 
gluon scattering amplitude should be taken with the extra factor 1/2 !. 

Accounting for the modification described we get the following result for four-fermion 
operators generated by graphs of the type given by Fig. 5: 

T2 = -^02, 02=(^EW5?) • (47) 

Dotted boxes in Fig. 4 present are reduced to local operators for fermionic momenta much 
less than k. We also account for the graphs which gives logarithmical dressing of the boxes, 



expli J (fxCint i ^ j I exp^-J d'^xd'^y 7]a{x)Dapix - y)rii3{y)^ 



(45) 

V-,ix)=0 
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i.e. {a ln(A;^/ Q^))" corrections. These corrections are due to the range of momenta between k 
and Q and diagrammatically are graphs presenting gluon exchanges between legs of effective 
four-fermionic operators as well as penguin-type fermionic loops for these operators. They 
were absent in two-loop considerations because of vanishing anomalous dimension of the 
relevant operator Opj (see Eq. (0)). Operators O12 have nonzero anomalous dimension. 

To sum up logarithmic corrections for the matrix element (7*|T|7*) we apply the standard 
machinery of the renormalization group. The operator basis of the problem consists of 
operators Oi, O2, Opj, Op2 (see Eqs. (|4^ , ^ H)). Their coefficients ci, C2, cpj, cp^ 
are functions of /i^ where fi is the normalization point. 

The set of renormalization group equations in one-loop approximation looks as 



_d_ 
d 



Nf{Ql) 



127r2 
1 



Cf 



3 ' 



d^' 



CFj 



127r2 



(Cl + C2 



/i 



d 

diJ? 

d 
dfi'^ 



Cl 



C2 



ai (2Nf + 1 



11 



TT 

3a;i 



-C2 



■Cl 



The solution for ci, C2 has the form 

1 



Ci(/i 



(ll/47r262^2) 

11 



+ 



11 



C2(/io) 



C2(/U^ 



1 



ci(/^o) 



11 



67r6o7i 
3 



C2(/Xo 



6vr6o7i 

ai(/ig 



ai(;U 



72 ■ 



1 + (ll/47r262^2) 
3 



27r6o7i 

where anomalous dimensions 71,2 are 



27r6o7i 

Cl(/io) +C2(/io 



Cl(/io) 



+ 



11 



C2(/io) 





""i(/^o)' 


72 1 




_Q;i(/i2)_ 





(48) 



71 



+ 



(49) 



Qi(/^o) 
ai(yu2) 



71 



+ 



(50) 



7l,2 




(51) 



Coefficients c^j and cp'^ are obtained then by a simple integration. Initial conditions for 
Ci(/i2) are set by Eqs. (|||, |^ at /x^ = A;^, 



ci(A;^ 



4 9\k^) 
3 A;4 



C2(A;' 



CFj{k'^) = cp^ik^) = . 



(52) 
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We need to calculate the value of Ci72(/i^ = Q"^) because with logarithmic accuracy the 
matrix element (7*|T|7*) is given by 



The result for Cp2[fi^ = Q"^) is rather lengthy expression 

1 



g2(e«e(2))-(ge«)(ge(2)) 



(53) 



9^/(^9) ^4 ^2(^2) ^252 ^ (11/4^2) 



+ 



15 



33 



3 27r6o7i ^n^l^ff I 1 + 



1-K- 



+ 3 



15 



+ 



11 



where 



27r6o7i Stt^I'jI / 1 + 72 



1 - K 



2 + 71 
^(1 



K 



(2+71) 



l + hoai{Q^)\n — 



aim 

The polarization operator 11 is given then by integration over k, 



n 



IGtt^ Jq2 

Coefficients of expansion in ai are defined by integrals of the type 



dk"^ 
Q2 'k^ 



-&oai(Q') 



(54) 
(55) 

(56) 
(57) 



where 5i arises from account of the second coefficient in the {3 function hi (see, e.g., ||2 
and equals to 

^1 = 



hi 



(58) 



The final formula for the expansion of 11 in powers of ai(Q^) has the form: 



n 



1447r2 7r26g + (11/471^) ^ 



X 



15 



+ 



33 



3 27r6o7i ^i^^hl-il) I + 



15 



+ 



11 



1 



27r6o7i Svr^fog^^ j i 



T(n + 2 + 7i + (5i 

r(3 + 7i + 5i) 
r(n + 2 + 72 + 5i) 



r(n+l + 5i) 



r(n + l + 5i) 



+ 



(59) 



r(3 + 72 + 5i) 

Assuming for the moment 7i = 72 = = 0, i.e. omitting effects of anomalous dimensions 
and two-loop /3 function, we get for coefficients a„ defined by Eq. (|^) 

n — >oo 13 



367r262 
12 



-M"(^ + l)! 



(60) 



Comparing this result with the Eq. (|3^) we see that the contribution of the four-fermion 
operators, or of the three-loop graphs, dominates at large n over that of the two-loop graph, 
it contains an extra factor n. Moreover, this conclusion is not modified if nonzero 7^ and 
5i are accounted for. Indeed, the dependence on 5i is universal for two-loop and three-loop 
contributions and drops off from the ratio of a„. As for the anomalous dimensions 7, of 
the four-fermionic operators (see Eq. (|5lD ), 71 is negative, 72 is positive. Notice that in the 
realistic case of QCD one also gets quite a few different diagonal operators both with positive 
and negative anomalous dimensions [^, ^ . In the large n limit the operator with the most 
positive 7 dominates so that account for the anomalous dimensions only strengthens the 
conclusion on the dominance of three-loop graphs. 

In the particular case when the number of flavors Nf is large, Nf ^ 1, one arrives at 
7i = —2, 72 = 0, 5i = and Eq. (|59|) leads to 

^ (-6or(r^ + l)! . (61) 



Although in large Nf limit the three- loop contribution is suppressed by an extra factor 1/A^/ 
as compared with the two-loop one it has an extra factor n. 

As an example of advantage of OPE technique used above let us mentioned that the sum- 
mation automatically accounts for rather complicated graphs, in particular, those where the 
high momentum gluon propagator is inside of fermionic bubble for a low (~ Q) momentum 
gluon. 

To summarize, we have demonstrated in this section that three-loop graphs generate 
new, four-fermionic operators. The technique of the operator product expansion works the 
same simple as in the two-loop case. Although the final expressions accounting for the 
anomalous dimensions are somewhat cumbersome they are straightforward. The most non- 
trivial conclusion is that four-fermionic operators dominate at large n. 



5 Calculability of ultraviolet renormalon. Conclusions. 

In this section we comment on the definition of the ultraviolet renormalon. The basic obser- 
vation is that if one tries to identify the ultraviolet renormalon as a set of graphs producing 
a certain asymptotic behavior, a„ ~ n\{hQ)^, then this set of graphs is in fact ill-defined. 

The demonstration of this point is straightforward. Let us go back to our operator 
product expansion (|T^. The coefficient functions are represented as series in aiik"^) 

Ci{k'^) = ho + hiaiik'^) + h2al{P) + ... + hia[{k'^) + ... (62) 



13 



The standard and tacit assumption is that one can approximate Ci{k'^) by, say, first term 
in this expansion as far as ailk"^) is small. This logic does not work however if we use 
the operator product expansion ( p!7D to evaluate the asymptotic of the coefficients a„ of 
expansion in our "original" as(Q^). Indeed, from Eq. ( |^ we conclude that for any finite I 
the contribution of hi to the asymptotic of a„ has the same n dependence. This means, in 
turn, that all the terms in the expansion (|6^) are the same important. 

Technically this happens because originally there are two large parameters, n and ln{k'^ / Q^) , 
and only the logs are controlled by the renormalization group. As a result of integration 
over k"^ powers of logs are converted into powers of n (see, e.g., Eq. (^I])) and the two 
parameters merge. The reason why all hi above contribute the same order of magnitude 
to asymptotical a„ is that lower powers of logs have larger statistical weights. This can be 
traced directly by inspecting combinatorics of I renormalon chains of graphs. As a result the 
two large parameters, combinatorial n and Infc^, get mixed up and the actual asymptotic 
of a„ is determined by an extremum in the two-dimensional parameter space. Finding this 
extremum looks much harder a problem than fixing the leading logs. It is not ruled out, in 
particular, that the contributions of various hi in Eq. ( |5^ cancel between themselves and 
the true asymptotic is very different from (^, (^. 

It is worth emphasizing that this difficulty with obtaining the asymptotic of the expansion 
coefficients is specific for the ultraviolet renormalon and does not plague the calculation of 
the infrared renormalon. The reason for this looks pure technical. On the other hand, 
renormalons are refiections of the Landau poles in the running couplings (see, e.g., |@]) and if 
indeed the true asymptotic of the UV renormalon is different from (^,(|^) then the Landau 
pole in ultraviolet does not manifest itself in, say, polarization operator n((5^). Thus there 
is a remarkable difference in the statuses of the infrared and ultraviolet renormalons. 

As another illustration of the uncertainties in the evaluation of the ultraviolet region in 
perturbation theory let us consider a tower of renormalons which means that we introduce 
another scale, k' such that 

{k'y > > Q\ {k'f ~ k'^e' (63) 
and / is large enough. Then coefficients hi are given by 

hi ~ /!(-M' • (64) 
The corresponding contribution of each hi to the asymptotic of a„ is again 

On — const ■ nl^—hoY . (65) 
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Moreover, we ignored the anomalous dimension of the operators which arise at the scale k' 
so that the true asymptotic can again be enhanced. But even these simple estimates suffice 
to demonstrate that no matter how large the virtual momenta are they do not detach from 
the asymptotic of a„. 

Thus, we conclude that at present there are no technical means to evaluate the true 
asymptotic of the UV renormalon. 

As for possible phenomenological implications of the results obtained we have unfortu- 
nately little to say. Pessimistically, one may argue that all the problems with evaluating the 
UV asymptotic are not physical and the reason why we fail is simply because we want to 
represent as a series in as{Q'^) the coupling ^^(fc^) which is much smaller than ^^(Q^) at 

^ Q^. To this end we need to keep many sign alternating terms absolute value of which 
is much larger than Q;s(fc^) itself. If we do not endeavor the expansion, however, the whole 
UV contribution is small and uninteresting. 

On the optimistic side, one might say that the divergences of the perturbative expansions 
reveal power like corrections and the failure to evaluate the UV renormalon signifies that the 
corresponding l/Q"^ terms are important. Indeed, applying the Borel summation one can 
convert the UV renormalon series into l/Q"^ term: 



N large ^ 



and if the constant K happens to be large numerically then the correction is also large. 
Moreover, one can generate the terms without invoking the Borel summation either by 
varying the normalization point of in the expansion or by exploiting the conformal 
mapping of perturbative expansions 0. 

One general remark concerning the power like corrections is now in order. Generally 
speaking, one would not expect these corrections be important numerically since they are 
naturally screened by lower order perturbative terms. The examples of successful phe- 
nomenological applications are known so far for IR renormalons only. For example, the 
ffist IR renormalon is associated with term 123]. Phenomenologically, the correspond- 
ing correction is described by the so called gluon condensate [^, (0|a;s(G'J^j,)^|0) where GJ^^, 
is the gluonic field strength tensor. The basic assumption is that this condensate is large 
and exceeds at moderate even the lowest terms of the perturbative expansion. Having 
only perturbative expansions in hand, this assumption is very difficult to justify. 

Similarly, to have the UV renormalon important we need to assume that the constant K 
is large. Better to say, we need a mismatch between large K in the asymptotic and numerical 
values of the lowest a^. The complicated structure of the UV renormalon revealed above 
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may be considered as suggesting such a complicated structure of the perturbative series. 
The main problem is not that such an assumption would be too bizarre but rather lack 
of a framework which would allow to develop a phenomenology basing on this assumption. 
Indeed, in case of the terms the main point was the possibility to relate these terms in 
various channels, that is not possible at present in case of the UV renormalon. 

Although we do not have constructive ways to introduce this phenomenology let us note 
that the natural direction would be Nambu-Jona-Lasinio type models since in these 
models one exploits four-fermion interaction as a low-energy effective interaction in QCD. We 
have seen that UV renormalon is naturally related to four-fermion interaction as well. One 
could hope to match 1/Q^ corrections due to the renormalon with NJL type phenomenology. 
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Figure captions. 

Fig 1. Building up the simplest renormalon-type graph. Dashed line denotes gluon while 
solid lines refer to fermions. One starts with an exchange of a vector particle of momentum 
k and inserts vacuum polarization bubbles n times. 

Fig 2. The graph used to evaluate coefficients of operator expansion associated with the 
simplest renormalon graphs of Fig. 1. Momentum k carried by the gluonic line is considered 
to be large. The fermion is understood to propagate in external electromagnetic and glu- 
onic fields so that the graph is in fact a subgraph of the one- gluon exchange depicted in Fig. 1. 

Fig 3. One-loop graph describing transition of electromagnetic current into a virtual 
photon with momentum q. 

Fig 4. Three-loop skeleton graphs giving rise to four-fermionic operators. Momentum 
k fiowing through the gluonic lines is considered to be large so that the operator product 
expansion is an expansion in inverse powers of k^'^. The dotted boxes mark subgraphs pro- 
ducing four fermionic operators O2 and Oi (see Eqs. (^Tf) and (H), respectively). 

Fig 5. The graph giving rise to the operator O2 in the limit of large k. It is a subgraph 
of the first graph in Fig. 4 where it is marked by a dotted box. 
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Fig. 3 



